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Abstract 

In reference ]|, a kinetic equation for gas flows was proposed that leads to a set of four macro- 
scopic conservation equations, rather than the traditional set of three equations. The additional 
equation arises due to local spatial random molecular behavior, which has been described as a 
volume or mass diffusion process. In this present paper, we describe a procedure to construct a 
Gibbs-type equation and a second-law associated with these kinetic and continuum models. We also 
point out the close link between the kinetic equation in [l|] and that proposed previously by Klimon- 
tovich, and we discuss some of their compatibilities with classical mechanical principles. Finally, 
a dimensional analysis highlights the nature of volume/mass diffusion: it is a non-conventional 
diffusive process, with some similarities to the 'ghost effect', which cannot be obtained from a fluid 
mechanical derivation that neglects non-local-equilibrium structures, as the conventional Navier- 
Stokes-Fourier model does. 
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I. INTRODUCTION 



In the continuity equations of fluid mechanics, as well as in the Boltzmann kinetic equa- 
tion for gases, spatial diffusive or source terms are not included. This is widely accepted as 
representing local conservation of mass. In fact, the continuity equation has a more general 
and purely mathematical foundation, originating from conservation of a probability (actu- 
ally, a probability density) that is substituted for rationally-defined physical quantities in 
order to derive a variety of well-known equations in physics. In conventional kinetic theory 
(and also in fluid dynamics), the spatial molecular probability density is taken to be equiva- 
lent to the physical mass-density of the fluid. Consequently, the mathematical conservation 
of probability density, i.e. the continuity equation, is taken to be the same as the evolu- 
tion equation for the physical mass-density. Mass-density at the macroscopic scale and in 
thermodynamics is taken to be an amount of mass divided by the volume containing that 
mass. 

Questions regarding the role of microscopic molecular spatial fluctuations and molec- 



ular/continuum length scale differences have often been raised by researchers 



5] . Some 



went on to suggest modifications to both the Boltzmann kinetic equation and the continuum 
fluid equation set, including dissipative components appearing in the mass-density equation 
y, y, 0-91. However, these models are often accused of violating the fundamental concept 



of mass conservation, and so they usually do not attract much interest within the fluid 
dynamics community. 

In rarefied gas dynamics, there are flow systems where experimental evidence indicates 



10 



151]. A number of these flow 



that existing gas flow models are not adequate, e.g. 
configurations (such as shock wave structure, sound wave propagation, the mass-flowrate 
in a microchannel) are therefore useful testing-grounds for any newly-developed gas flow 



models. Sone et al. 



16 



19| described a series of flow configurations that, although belonging 



to the standard Navier-Stokes-Fourier regime, cannot be described by the conventional fluid 
model. With increasing interest in micro/nano scale flows, experimental evidence from 
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both liquids and gases [20|, |2l| is extending the range of apphcation domain in which the 
conventional fluid model fails. There is, therefore, a substantial motivation to develop and 
test new flow models in order to understand the physics of flows beyond conventional fluid 
mechanics. 

In the kinetic theory of gases, there are also recurrent problems faced by continuum mod- 
els derived using approximation solution techniques to the conventional Boltzmann kinetic 
equation to an order in Knudsen number beyond that of Navier-Stokes-Fourier. These are 



22|: consistency with the second law of thermodynamics, frame indifference, and angular 



momentum conservation. Woods extensively investigated these problems for Burnett-type 



equations 



23-25 



In reference 



l|, a kinetic model equation was proposed, followin g a claim of the existence 
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-|28j. In this present 



of a new form of diffusive process, so-called "volume diffusion" 
paper, we show that it is possible to construct a Gibbs-type equation and a second-law 
associated with that kinetic model. Then, using a dimensional analysis, we demonstrate that 
volume/mass diffusion terms in the fluid equations are non-conventional diffusive processes 
reflecting the microstructure behaviour beyond the accuracy order of Navier-Stokes-Fourier. 



II. STOCHASTIC KINETIC EQUATIONS AND CONTINUUM-FLUID VOLUME 
DIFFUSION 

A. The Klimontovich spatial stochastic kinetic model 

The Boltzmann kinetic theory, based on a single particle distribution function, attaches 
no special importance to the spatial arrangements of molecules. An average number of 
molecules is associated with a position variable X. This position variable is then associated 
with the continuum fluid position variable. No clear distinction exists between the kinetic- 
level spatial position and the hydrodynamic position variable. Klimontovich questioned this 
aspect of the kinetic description, and introduced a separation scale in these two position 



variables 



7|. He assumed a microscopic distribution of the continuum volume given by 
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H3l|. 



a Gaussian. Then, using some relaxation arguments and integration over the microscopic 
volume, arrived at a modified Boltzmann equation which differs from the original by a spatial 
diffusion term. Klimontovich's modified kinetic equation can be written: 

df 

^ + e ■ V/ - V ■ «:V/ + ■ - I^W) = 0, (1) 

where k is a spatial diffusivity coefficient, ^ is the kinetic particle drift velocity, F^xt denotes 
an external body force (not included originally in the Klimontovich expression), is the 
gradient operator in velocity space, I^{f) is the Boltzmann collision integral, with f{t,X,C,) 
the (modified) distribution function of the molecules. In fact, equation ([1]) is a Fokker-Planck 
type of equation for diffusion with drifts, where the diffusion is explicit in both position and 
velocity spaces, and the velocity space diffusion operator is substituted for the Boltzmann 
collision integral. The equation can also be derived using purely stochastic reasoning 
On dropping the spatial diffusion term, it reverts to the classical deterministic Boltzmann 
equation. 

A macroscopic continuum set of equations corresponding to kinetic equation ([T]) can be 
derived by first defining macroscopic fiow variables as follows. The macroscopic mass-density 
p{t, X) and macroscopic fiow velocity U {t, X) are 

p{t,X) = j Mf{t,X,Od^, and p{t,X)U{t,X) = j M^f{t,X,Od^, (2) 

where M is the molar mass. Note that in equations the integration is only over the 
microscopic molecular velocity ^. The position variable X is the same as that appearing in 
kinetic equation ([1]). But, in the Klimontovich modification, variable X in both equations 
([1]) and ([2]) represents a continuum position variable, not a position at the kinetic level of 
individual molecules. 

The derivation of a continuum set of local transport equations follows by multiplying the 
stochastic kinetic equation ([T]) by M, M^,M^^/2, and integrating over velocity space only. 
This gives, respectively, 

^ + V ■ [pU -nVp ] = 0, (3a) 



dpU 



V ■ [puu] + V ■ [pi + n] -V ■ [kV {pu)] = 0, 



(3b) 



d_ 
di 



-pU^ + pe. 



+ V 



V ■ [{pi + n) ■ ?7] + V ■ [q] 



-V- 



0, 



(3c) 



where I is the identity matrix and p is the kinetic pressure, which is related to internal 
energy by 3p = 2pein- The underlined terms in this set of equations are the new additions 
arising from the random change in molecular positions at the microscopic level, not from 
the random exchange of momentum. Quantities e^, 11 and q are all related to the random 
component of molecular velocity in the same way as in the traditional derivation from the 
deterministic Boltzmann equation. Essentially, the shear stress 11 and heat flux q that result 
from molecular level diffusion of momentum and energy can be given their Newton and 
Fourier law expressions. According to the kinetic theory definition, the kinetic temperature 
T can also be associated with ej„ = (3/2) i?T, with R the specific gas constant. 

The set of equations ( l3a|) to (!3c|) produces, for example, a better prediction of sound wave 
dispersion, which is a flow configuration where microscopic structures are not as negligible 
as the Navier-Stokes- Fourier model requires 32| . However, the same underlined terms above 



that generate acceptable flow results, also appear apparently to contradict classical mechanic 
conservation principles and the second law of thermodynamics (see Appendix IVIip . 

Arguments for the modified kinetic equation ([1]) focus on the introduction of the spatial 
distributions of molecules. Assuming a Gaussian distribution of the continuum volume, 
as Klimontovich did, avoids introducing an evolution equation for the continuum volume. 
However, a microscopic variable could be introduced to describe the actual evolution of the 
molecular spatial distribution. Such a full derivation is presented in jlj], and we now analyze 
the thermodynamic consistency of that volume model. 
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III. EXPLICIT VOLUME-BASED KINETIC AND CONTINUUM EQUATIONS 



A. The kinetic and continuum equations, and mechanical consistency 

We follow a kinetic description with a distribution function f{t, X, ^, v) that describes the 
probability that a molecule at a given time t is located in the vicinity of position X, has 
its velocity in the vicinity of ^, and has around it a microscopic empty space given by the 
additional variable v. The presence of the microscopic variable v is justified by the stochastic 
aspect in the microscopic picture of the spatial distribution of molecules. It is chosen here 
as a microscopic volume in order to track the macroscopic volume occupied by the entire 
collection of molecules. 

Then, admitting that the probability density f{t, X, C,, v) must be conserved, we can de- 



rived a kinetic equation for f(t, X, ^, v) as 



H + (e ■ V)f + £ + F,., ■ V^f - /e(f) = 0. (4) 

Operator /^(f), describing the molecular level diffusion of momentum, is taken in equation (jl]) 
as the Boltzmann hard sphere molecular collision integral; this is because it operates only in 
velocity space, while X and ^ are independent in the same way as in kinetic equation ([1]). On 
the left-hand-side of equation (jl]) the term containing W now replaces the Laplacian spatial 
diffusion term in kinetic equation ([T]). This term embodies the evolution of the macroscopic 
continuum volume from the microscopic volumes, and W is the local continuum volume 
production rate. An expression for W needs to be defined, and it is a property belonging to 
a collection of molecules rather than an individual molecule. The external force Fext, which 
operates as a drift in velocity space and generates momentum and accelerates a molecular 
motion, will be neglected for simplicity without changing any aspect of the fundamental 
description. 
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1. The macroscopic field variables 



The introduction of an additional microscopic parameter modifies the method of con- 
structing macroscopic field variables in our explicit volume-based kinetic description. 

As f (t, X, ^, v) is presented primarily as a probability density function, we have a normal- 
ization factor (or a reduced probability density): 

A„(t,X) = j J f{t,X,tv)d,d^ . (5) 

The mean value, Q(t,X), of a gas property Q is then defined by, 

^(^'^) = ^TT^^ / [ QKt,X,^,v)dJ^ . (6) 



We have then a local average of v as the mean empty volume v{t,X) around each gaseous 
molecule, i.e. 



= Ajt;x) j j ""^^^^ ^' ^' ''^'^^'^^ ■ 

From equation ([5]) it appears that An represents a probable amount of mass (or, equiv- 
alently, number of molecules) around position X, but does not contain information about 
the distribution of that mass around that position. How the amount of mass is spread out 
in space is given through v{t, X) in equation ([7]). As average quantities defined via equation 
([6]) are mass-based averaged, a definition of the physical mass-density of a fluid should be: 

*-^' = ^- 

where M is the molecular mass. 

Two mean velocities are defined using two different weighting values: the local mean 
mass-velocity Um(t,X) is given through: 

An{t,X)UUt,X) = j j ^f{t,X,^,v)d^d,, (9) 

and a local mean volume- velocity Uy{t,X) by using the microscopic empty volume as the 
weighting, i.e. 

v{t,X)Ar,{t,X)U,{t,X) = j jv^f{t,X,C,v)d^d,. (10) 



It is important to note that the two definitions for and Um coincide if v is assumed 
constant. For a homogeneous v the present description is therefore equivalent to the classical 
kinetic theory. For a non- homogeneous v, we denote U-u — Um = v~^3i, which evidently has 
dimensions of velocity. The quantity Ui, — Um can physically represent either a volume or 
mass diffusion process 



2. Consistency with mechanics 



The separation of mass from volume in this kinetic construction introduces two types 
of averaging that relate to the way balance equations are analysed in classical continuum 
mechanics; namely, 'control volume' and 'control mass' concepts 33|. So the macroscopic 
set of derived continuum equations can be written in two versions. 

By integrating the kinetic equation over and v, after multiplication by the different 
microscopic elements v, 1,^ and the set of macroscopic equations take the form: 

dAnV 



dt 



+ V ■ [AnVUm] + V ■ [AnJy] = A^W, 



(11a) 



dA. 



dt 



- + V ■ [AnUm] = 0, 



;iib) 



dA^U„ 
dt 



+ V ■ [A^UmUm] + V ■ = 0, 



:iic) 



d_ 
di 



-AnU^ 



d_ 
di 



[Ant-.n] + V 



-A^U^Um ~\~ A-i^tinU^ 



+V ■ [A„^ ■ f/„^] + V ■ [A„.q] = 0, 



;iid) 



where Um is the average mass-velocity, and quantities J^, etn and q are all quantities 
remaining after extracting their macroscopic parts due to Um] so these quantities are all 
related to the random mass velocity ^ — Um- At this stage, the set of equations (lllbp to 
fllldl) has a strong conservative form related to mass; in fact, they are exactly the same as 
the traditional expressions for local conservation of mass, momentum and energy. So there is 
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no violation of local angular momentum conservation, or any other fundamental mechanical 
problem. Separately, the actual volume of fluid containing the local mass, momentum and 
energy is described in a separate equation filial) , which does not have the conservation 
structure of the others because of the production term W. 

The volume of a fluid is generally considered as a thermodynamic property, so its variation 
should impact on other properties too. The new flux of volume, ^"^J^, has the dimensions 
of velocity. Adding this velocity to the average mass motion Um in order to define convective 
fluxes, we obtain a second form of the momentum and energy equations that accounts for 
the fluid volume variations. The macroscopic set of equations filial) to fllld|) is therefore 
equivalent to: 

DA„ 



Dt 



(12a) 



A, 



Dv 
'Dt 



-V ■ KJ,] + A„iy, 



(12b) 



A, 



DUrr 



-V ■ A J P' 



1-^ 



(12c) 



2^m I "--jn 2v'^ 



V-An 
V-Ar, 



p' T T 



(12d) 



where the material derivative is defined with respect to Um as D/Dt = d/dt + Um ■ V, and 
Uv = Um + v~^3y. In equations fll2ap to fll2dp the fluid internal properties and fluxes ac- 
counting for not only macroscopic mass motion but also the variation of its volume occupied 
become P' and q' for momentum and energy flux, and e-„ for the internal energy. It is also 
clear that the structure of the local transport equations for momentum and energy is now 
different. Taking P' and q' as diffusive fluxes, they should follow the classical phenomeno- 
logical Pick's law, at least to a first order linear approximation. These are derived later by 
considering compatibility with the second law. 
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Taking equation fll2cp to the cross product with a continuum position vector X we write 



-X A V ■ A„ ( P' - —J, J, , , 



(13) 



which is equivalent to 

An^^ [X A f/„] = -X A V ■ An (P' - ^ J. J.^ . 

Now, for a symmetrical second order tensor T, the following property holds: 



(14) 



X A 



V-T 



X AT 



(15) 



Pressure tensor P' is symmetric by definition (see its full expression in equation I28ap . Con- 
sequently, sum of tensors appearing on the right-hand-side of equation (fT4|) is symmetrical. 
Equation f|T^ has therefore the following final form 



An^ [X A U„ 



AnX A P' - — J,J, 



(16) 



Equation f ll6p shows that the volume-based momentum equation fll2cp still satisfy the prin- 
ciple of angular momentum conservation. Generally, because the velocity vector involved in 
continuity equation f ll2ap is the same as that involved on the left-hand-side of momentum 
equation f ll2cp . the set of volumed-based equations fll2ap to (112dp satisfies following me- 
chanical properties: galilean invariance - integrability - angular momentum - center-of-mass 
position; see explicit demonstrations in Appendix IVIIII 



B. Consistency with the second law 



Evaluation of the second law of thermodynamics in this case requires an energy equation 
that accounts for volume variations. First, the energy equation (I12dp is transformed using 
the momentum equation fll2cp into: 



D 
Dt 



V-Ar. 



e -1jM-J 

I V 



(17) 









V -An 


q + P' ■ - J. 

V 





/1 .2 
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Then, using the continuity equation f ll2ap and pP' = p'l + IIv, with p' the kinetic theory 
definition of pressure by, 3p' = 2pe^„, we obtain 



A, 



D 



p' —1^ 



-V- A 



p Dt 



q + P' ■ - J. 

V 



+ V- A„ 



p' -— T , . 



(18) 



p 



The volume evolution equation (112bp may be rewritten in terms of the density defined in 
equation as. 



ID An DAnp-' ,1 ^ r, T 1 



which, substituted in equation flTSl) . leads to: 



(19) 



™ 2^2 



P 



Dt 



■■^V.[AnJ.]-p'^W 



V- A 



-W-An 



I V 



q' + P' ^Jy 

V 



An { - -^J,J, 1 : Vf/„ 



(20) 



This can be rearranged using = 3/2(p'/p) into : 



A, 



D 

'Dt 



e'. 



2v 



,DAnP 



p 



Dt 



p 



'^W - V ■ 

M 



A 

pv 



Hv J, 



+ ^J,J, :V?7„ (21) 



P 



+ 



f J V 



V- 



2 p 

—'n.Um '■ — Hjv : Vf/m , 

P P 



where the final three terms arise because the shear stress has been split into its components 
due to Um and due to J^, respectively. Hum and 11 jv, such that IIv = ^Um + IIjv, as 
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From the structure of the energy equation fl21l) . we define an entropy quantity s such 
that: 



A, 



D 
Dt 



,DAnP 



p 



Ar,T 



(22) 



Dt Dt' 

with T' the associated temperature. The entropy equation obtained after substituting this 
definition into the energy equation fl2ip is, 



Ds 

m 



-p 



M 

p' 



pv 



-Hv J, 



+ 



^ J J _ 

p 

AnP' , An -J 2 



Jv 



: Vf/„ 



(23) 



P 



+ ^2 "^^ 1 



V- 



■ 3Anp'l' 
^nQ + 77 — ~ — ~'Jv 

2 p V 



A 



P 



We have here a non-local-equihbrium entropy defined via equation fl22i) . which is a Gibbs-hke 
equation, and an entropy evolution equation (l23l) in which the volume production rate W is 
assumed to be a function of the macroscopic thermodynamic properties, but with (so far) an 
unknown constitutive equation. It follows that the first and second laws of thermodynamics 
can be used to construct an equation for W from equation (!23|) . 
First, note that if we set 



M 















V- 


— n J 


+ 




: Vt/„ + — Hjv : V 






pv 




^V^ J 


P 


V 



(24) 



+ —Hum : V 
P 



V ■ [A„J^] + V 



^nP' , An J 2 



—9 



the entropy equation fl23p becomes 



AnT 



Dt 



or equivalently. 



-V 



2 p f 



^Hf;^ : VU,n - — Hj, : V 
P P 



(25) 



AnT' 

Dt 



■ ^AnP'l' 
-^n^. ~r — - —Jv 

2 p V 



— Hv : VUv. 
P 



Equation 026 p has a structure such that entropic heat fiux qg is given by, 

_ r , 3 Anp' 1 ■ 

L 2 p V 



(26) 



(27) 
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According to Linear Irreversible Thermodynamics concepts, and following equation fl26p . 
the first order linear representation of the diffusive shear stress and entropic heat flux can 
be given according to a Pick's Law type of representation. These are then written, 



ij. = -^Vv = ^Vp , (28c) 

V V p 

where = Um + v~^3v The coefficients /i', k'^, t]' and respectively, dynamic 

viscosity, heat conductivity, bulk viscosity, and the mass-density diffusion coefficient in this 
volume-based description. As the kinetic pressure p' was defined by the trace of the pressure 
tensor, we also have (2/3)/i' — 77' = assuring the symmetric-traceless property of the shear 
stress tensor. The entropy evolution equation fl2B]) may now be rewritten: 



A T'— - V ■ 



- yHv : VU^. (29) 



P 

In summary, with an entropy evolution equation f l29p and the volume production term given 
by equation ([21]), we have ensured that the set of volume-based continuum equations is 
thermodynamically consistent with the first order linear expressions for shear stress and 
heat flux. Indeed, the structure of the entropic heat flux in equation fl27|) is such that the 
transport coefficients satisfy the Onsager reciprocal relations, and the entropy production 
rate follows a bilinear structure. Then non-negativity of the coefficients ^ A^j p and t^'^^An/ p 
assures the non-negativity of the entropy production with a shear stress expressed as in 
equation (I28ap . However, we note that, to satisfy these properties, constructing W as it is 
given in equation flM|) is not a unique option. 

The expression for the volume production, equation (l24l) . can now be substituted into 
the volume equation f ll2bp to finalize the evolution equation for the specific volume of the 
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fluid: 



A, 



p' Dv 













— n J 


+ (^J.J.) : Vf/„ + y : V 






pv 




V 



(30) 



V 



AnP' , A 



p 



l2 "^v 



It is easily seen in this equation that a vanishing volume flux also means a constant v in 
the material reference frame. In such circumstances, all additional terms due to volume 
evolution in the momentum equation fll2cl) . as well as in constitutive equations f l28al) to 
5c]), cancel out. The set of volume-based equations then simply returns to the traditional 



Navier-Stokes- Fourier set, with a mass-density represented by MA^, 



IV. THE NATURE OF THE NEW DIFFUSIVE VOLUME/MASS TERMS 



A. The problem of a fluid in rotating equilibrium 



ation of mechanical properties by 



35| . The conflguration is a fluid or 



A conflguration which was used to point out the vio_ 
volume/mass-diffusion hydrodynamic models is given in 
gas in rotating equilibrium, at constant temperature. More precisely, for a fluid rotating in 
equilibrium at constant temperature, contributions from VT' and the symmetric-traceless 
Navier-Stokes shear stress component related to Vf/m vanish. Entropy equation (129]) . con- 
sidered without any preliminary analysis of the nature of the terms involved, can then be 
reduced to: 



, ^,Ds An„ „ 
Dt p 



V 



(31) 



Equation flM]) shows a production of entropy for an apparently equilibrium conflguration. 



This is, according to 



35|, not acceptable. 



In fact, whenever additional local microscopic structures are involved, additional en- 



tropy production in the way it appears in equation fl3T|) s 



methods for non-local-equilibrium fluids, for example 



36 



lould be expected (see existing 



39|). We now demonstrate us- 



ing a dimensional analysis that the production term in equation ( ISTj) is indeed a type of 
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non-local-equilibrium term which is negligible at the order of approximation of the Navier- 
Stokes-Fourier model of hydrodynamics. 



B. A dimensional analysis of the new entropy (energy) equation 



Denoting the gas mean free path A, we define the following reference fluid parameters: 



,,0 _ J^r^x ^.0 _ - 

fi — p U A, l^m — ~y ^ — 



(32) 



,0 _ 



P° = MAl 



where superscript refers to a reference characteristic flow parameter, and C° is a charac- 
teristic molecular agitation speed. Then we define dimensionless flow properties through: 



* 
/i = /i 



* /lO 



.* .1 *0 /I /I* /I 



(33) 



where superscript '*' indicates a dimensionless quantity, and L a macroscopic length scale. 
Entropy equation fl2^ can then be rewritten 



Ds* K°r° 



L2 



-V- 



A 



(34) 



02^0 ^ 



p* 

p* 



p 



n r -* ■'-'■Jv ■ ^"^n 



which can be re-arranged using the expression for p° and from equation f l32|) as 

'A! 



P' 



0^02 to 



Dt* 



p* 



Um 



A* 

p 
p* 

^Vp* 
p* 



vu* 



p* 



(35) 



in which appears the Knudsen number, Kn = X/L. Defining the reference entropy per unit- 
mass s^, in a way that the unit- volume of fluid is re-scaled by the mean free volume A^, 
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I.e., 



f_ _ 1 

entropy equation (l35|l then yields the following full dimensionless form: 



(36) 



Ds* 



Dt* 



A* 

flu 

p* 



A* 

p* 



(37) 



A* 

p* 
A* 



Urn 



: V 



''-Vp* 



A* 
p* 



Jv 



p* 



■Vp* 



In equation fl37p it appears that all volume/mass diffusion terms are of higher order than 
Navier-Stokes-Fourier. The term related to the entropy production observed in the rotating 
fluid problem in equation ( 13T]) . is of order K^. 

In kinetic theory the Knudsen number is the parameter that quantifies the degree of non- 
equilibrium. It is considered as a small parameter when deriving continuum hydrodynamic 
models, as in the Chapman- Enskog expansion technique. So truncating equation ( ET^ at 
order K^, we get the zero entropy production that is consistent with the Euler hydrodynamic 
regime. Truncated at order Kn, equation ( 137|) becomes exactly the Navier-Stokes-Fourier 
entropy evolution equation. All volume/mass diffusion terms appear as non-equilibrium 
effects at an approximation level higher than that of Navier-Stokes-Fourier. 



Sone et al. 
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19[ have described a series of flows that, while belonging to the continuum 
fluid regime, the conventional fluid mechanical model is incapable of predicting. Instead, 
one has to rely on additional higher order Knudsen number terms - traditionally called 
rarefaction regime terms - in order to predict the flow behaviour. The dimensional analysis 
above presents a clear exposition of these intriguing situations. Indeed, in the derivation 
of the volume/mass diffusion model we did not place any restriction or assumption on the 
Knudsen number. So, there may be configurations where these high order terms may dom- 
inate the lower order terms. This may be the case for flows with A'„ > 1. In other words, 
our volume/mass diffusion terms are not of the same type as the conventional macroscopic 
Navier-Stokes-Fourier diffusion terms; they are effects arising from another level of micro- 
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scopic behavior that may affect the global ffow structure (or may not) depending on the 
configuration. In the case of the rotating fluid in equilibrium, they do not. Derivations of 
fluid mechanical models that neglect local microscopic structures, as many do cannot cap- 
ture these types of physical effects as they simply do not account for the material structure 
at that level. 



V. THE VOLUME/MASS DIFFUSION MODEL IN AN APPARENT CONTIN- 
UUM LIMIT REGIME 

From our analysis of the nature of diffusive volume/mass terms, it appears that to account 
for only the first level of this non-local-equilibrium embodied in the new diffusive volume 
process, we should naturally neglect terms of high order J^, i.e. assume terms such as 
3^3 V to be negligible. In so doing, we may also argue that at this first level, while a detailed 
evolution equation of the new local microscopic volume may not be important by itself, it still 
effects the macroscopic flow, mass, momentum and energy equations. In such circumstances. 



equations fll2ap to fll2dp may be reduced to a three set of equations written: 



Dp 
Dt 



-pV ■ U„ 



(38a) 



DUm 



(38b) 



-U^ + e' 



closed with 



-V ■ [p'U, + Hv ■ f/.] - V 



' / 3 ,1 
2 V 



P 



2 pv 



Uy = Um + -Jy = Um + — Vp, 

V p 



(38c) 



(39a) 
(39b) 
(39c) 



where MAn = p = M/v is identified with the fluid physical local thermodynamic density 
and the evolution equation for the microscopic volume disregarded. This set of equations is 
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exactly the same as proposed more recently also by Brenner 



28 



40|, using entirely classi- 



cal irreversible thermodynamics reasoning. The set of equation s fl38a 



on s fl38a 
|9|,fl34 



-fl39cl) . as opposed to 

other previous attempts of volume or mass diffusion models [9|, |32|, |3j, |4l| can be said fully 
mechanically consistent (see explicit demonstrations in Appendix IVIII|) . Note that, mass 
conservation equation fl38ap does not involve explicitly a diffusive mass flux component. In- 
stead, the volume/mass diffusion appears via the energy equation and in the expression of 
the shear stress tensor. If the volume velocity from equation fl39cp is substituted for the 
mass velocity in equation f l38ap . then we observe evidently a form of diffusive component 
in the mass equation. However, we note that this is not a true physical 'dissipation' or 
'disappearance' of the mass as that diffusion is regarded with respect to the volume ve- 
locity rather than the mass velocity. Our volume-based momentum equation (138bp , closed 



with equation f l39ap . is also the same as momentum equation derived recently in |39|. In 
that paper, the authors presented the derivation using a different explicit spatial stochastic 
approach but did not provide the complete accompanying form of their energy or entropy 
evolution equations. 

Finally, the two kinetic model equations ([T]) and introduce into the kinetic description 
a stochasticity in position that does not appear, at least explicitly, in the Boltzmann kinetic 
approach. Equation ([1]) incorporates this only as a diffusive volume, while equation (jl]) 
incorporates it in a more complete picture as a convected volume and a volume source term. 
Consequently, equation ([T]) may be regarded as a reduced form of equation (jlj) and both 
kinetic equations relate to a description in which local-equilibrium does not necessarily hold. 

Without the volume source term in equation fl23|) . it would be challenging to construct a 
second law analysis. As a matter of fact, constructing a second law for high order hydrody- 
namics beyond that of Navier- Stokes- Fourier while maintaining their mechanical consistency 



22, 



24 
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42| . It appears that 



(e.g., for the Burnett equations) has always been a challenge 
this new theory of volume diffusion and the distinction between different types of velocities 
allow to achieve this. The procedure described here to obtain a version of the second law 
for this volume-based approach (in which the additional volume/mass terms are of Burnett 
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order), therefore opens the door to reinvestigating longstanding problems with most, if not 
all, extended hydrodynamic models. 



VI. CONCLUSION 



In this paper we have successfully constructed a thermodynamic second law for a volume- 



based hydrodynamic model presented first in 



l|. We have also provided further analysis of 



the nature of the volume/mass diffusion process that is still currently a subject of active 
debate. We have demonstrated that this diffusive process is not a conventional diffusion of 
the Navier-Stokes-Fourier type; it arises from a deeper level of microscopic behaviour that 
may or may not affect the global flow structure, depending on the flow configuration. Gen- 
erally, a derivation that neglects local microstructure, as most conventional fluid mechanical 
derivations do, cannot capture this. 
Future work may include: 

• reinvestigating thermodynamic problems related to Burnett and other extended hy- 
drodynamic equations; 

• comparing and contrasting the various volume/mass diffusion models proposed so far; 

• investigating 'ghost effect' flow configurations, which are flows in the pure continuum 
regime that require high order Knudsen number terms to be predicted appropriately 
- these are therefore an ideal testing ground for volume/mass hydrodynamic models. 

VII. APPENDIX: THERMODYNAMIC AND MECHANICAL CONSISTENCIES 
OF THE KLIMONTOVICH VOLUME MODEL 

Let us assume a 'material' derivative as, D/Dt = d/dt + U ■ V. Then equation ( l3b|l can 
be rewritten, 

DU , (dp j^^A ^ ^ y . ^ n] _y . [^y (^f;)] ^ g, (40) 



Dt \dt 
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which becomes, after introducing the mass-density equation ([3a] 
DU 



Dt 



V ■ [pi + n] + (V ■ [kWp]) f/ - V ■ [kV (pf/)] = 0. 



(41) 



The underhned terms are again, the additional terms due to the adoption of the Khmontovich 
kinetic equation. Taking the cross product of equation fHTl) with a hydrodynamic position 
vector X, we notice quickly that the underhned terms generate. 



X A {(V ■ [kWp]) f/ - V ■ [kV (pf/)]} 



(42) 



which we cannot write in local conservative form, i.e. as V ■ [...]. So these terms appear 



34| for other families of 



as local angular momentum production terms, as pointed out in 
dissipative mass continuum models. Equation fHTj) is therefore said to violate local angular 
momentum conservation. 



To analyze the second law of thermodynamics, energy equation fl3cl) is first re-written: 



P 



D 
Dt 



+ V ■ [{pi + n) ■ f/] + V ■ [q; 



kV -pU^ + pCi 



0, (43) 



which becomes, after introducing the mass conservation equation fl3al) . 



D 



+ V ■ [nVp] -U^ + e 



V ■ [{pi + U)-U] + V-[q 
1 



-V ■ 



0. 



(44) 



Momentum equation fHTl) can be used to eliminate the macroscopic kinetic energy terms; 
introducing the density equation (l3a|) . the energy equation is finally 



Dcin Dp ^ p 
p—p^ + PP^^ — 2kVp ■ Vein - KpV ■ Vein + -V ■ [^VpJ 

-KpVU : VU + n : Vt/ + V ■ [q] = 0, 



Dt 



(45) 



in which ':' denotes the Frobenius inner product. In classical fluid dynamics, the specific 
entropy s is defined by adopting the Gibbs (local equilibrium) equation: 



Ds Dein , Dp-' 



(46) 
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Using equation (H^ . energy equation (HSj) becomes an equation for the entropy: 

P^-^V-V(pT) + «:(i? + OV-Vp 

-Kp^VU : Vf/ + : Vt/ + ^ V ■ [q] = 0, (47) 



where we have also used ej„ = c^T . Finally, with the identity, 

^v-v.^^.v.(^), («) 

where is a scalar field, the entropy equation (H7I) takes the form: 

Ds ^ /V(pT)\ fVp\ „ /VT 

P;^-P«^c,V ■ (^^^J + p/€ (i? + c„) V ■ (^-^J - ■ 

'pU-.VU - : + ^^VT . VTH- {^Vp . VT - ^Vp ■ Vp}. (49) 

From this entropy equation (H9|) we observe that the entropic heat fiux is — (k/j + pKCt,)VT, 
which shows that heat conductivity in this stochastic continuum model cannot be associated 
systematically with the classical heat conductivity. The terms in the curly brackets in 
equation fj^9|) can be either negative or positive. This suggests that a negative temperature 
or decreasing entropy can occur. According to classical thermodynamics with the Gibbs 
equation fj46l) . these terms are therefore undesirable. Entropy fiux is given on the left hand 
side of the equality in equation (jinj) by terms associated with the divergence operator. The 
unwanted and difficult terms in both the momentum and entropy equations are generated 
by the diffusive term in the density equation (or, more precisely, the spatial diffusive term 
in the initial kinetic equation). 

However, in making this conclusion the fact is ignored that, in this Klimontovich descrip- 
tion, the hydrodynamic position variable X represents the motion of not a given molecule 
with a precise mass but a collection of molecules occupying a volume that may be chang- 
ing. More precisely, the velocity U appearing in the set of equations fl3aj) to fl3cl) may be 
attributed to a volume velocity in same way as presented in the set of equations fll2ap 
to fll2dp . Then, because conservation properties such as mass conservation are all prop- 



erties related to the mass and mass velocity /current, we may conclude that the violation 
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of angular momentum conservation in equation, is simply an apparent inconsistency: this 
is because within equation fj4T]) and the derivative d/dt -\- U ■ V, we have been identifying 
erroneously the volume velocity with the mass- velocity leading to an incorrect interpretation 
of the violation of the conservation rule. A clear picture is given in the next appendix. 

VIII. APPENDIX: VERIFICATION OF OTHER MECHANICAL PROPERTIES 
OF THE FULL VOLUME-BASED SET OF EQUATIONS 

Mechanical and thermodynamic properties to be satisfied by a hydrodynamic model as 



described and listed in 



34| are: 



galilean invariance 



integrability 



angular momentum conservation 



steady rigid fluid rotation 



center-of-mass position 



second law of thermodynamics 



We take the set of equations f l38ap - fl39cp in which the material derivative is with respect 
to Urn- Note the conservative form of that set of equations as: 



1^ + V ■ [pU,n] = 



(50a) 



d_ 
di 



[pu^m] + V ■ [pUM + v-[pi + n,] = 0, 



(50b) 



di 



pU^ + pCin U, 



+ V ■ [[pi + Hv] ■ f/„] 



V ■ i-KhVT) = 0, 
(50c) 
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where for simplicity in notation we drop the prime above variables p', e-^, T' and also bar 
above p. Then, according to equation flSOap the mass flux is given by pUm and it is the 
same as the momentum density appearing the left-hand-side of the momentum transport 
equation fISObI) . 



A. Galilean invariance 



We consider the following transformation: 



t = t 

X-ct =X , (51) 

where c is a constant velocity vector. Equations flSTl) give transformed partial time and 
position derivatives by, 

a 



d_ 

at 
a 

ax 



at^^ ax 
a 

ax 



(52) 



which gives the transformed material derivative as: 

D d 



f/™ ■ V 



d_ 

dt 



D 



(53) 



Dt dt 

so the material derivative is galilean invariant. Then, substituting this material derivative 
and the change of variable (15T]) in the mass and momentum equations fl38ap and (]38bp closed 
with equation f l39ap . we obtain the same form of equations. The galilean invariance of the 
energy equation is seen in a similar way using its expression written in (!65|) and also due to 
the form of diffusive fluxes in equations fl39al) - fl39cp . and also as transport coefficients are 
generally assumed independent of position and time variables. 



B. Integrability 



Integrability is concerned with the transport equation for total mass flux. In the model 
of equations fl5Ual) - fl50cp the total mass flux is pUm- It follows that the transport equation 
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for the total mass is just the momentum transport equation (150bP : so it as a conservative 
equation. 



C. Angular momentum 

Conservation of angular momentum is one of the most important properties. Taking 
equation (138bl) to a cross product with a hydrodynamic position vector X: 



X Ap- 



Dt 



-X A V ■ [pi + Hv] 



(54) 



which is equivalent to 

p-^ [X A U,n] = -X A V ■ [pi + Hv] . (55) 
Now, for a symmetrical second order tensor T, the following property holds: 



X A 



V-T 



X AT 



(56) 



The pressure tensor IIv as given in equation f l39ap is symmetrical. So equation (155|) has the 
following final form: 



: -V ■ [X A [pl + Hv]] . 
Therefore conservation of angular momentum principle is satisfied. 



(57) 



D. Center-of-mass position 

In absence of external force the center-of-mass motion must be uniform. To this end, one 
expects the quantity. 



B = pX- pUmt 



(58) 



to be a conserved quantity. So we write 
dB d 



do d 

dt ^[p^-pu^t] = x^-t-[pu^]-pu„ 



(59) 
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in which transport equation of mass and momentum equations flSUal) and (150b|) can be 
introduced to obtain, 

dB 

— = -XV ■ [pUm] + tV ■ [pUmUm + [pl + H^]] " pU^. (60) 

Equation (plj) can be written, 
OB 

— = -V ■ [pXU^] + pUm + tV ■ [pU^Um + [pl + Hv]] - pUm, (61) 

or simply 

— = -V ■ [BUm -t\pl + Hv]] , (62) 

which is therefore a conservative transport equation for the quantity defined in equation 
( 158|) . Consequently center-of-mass position principle is satisfied. 

The fundamental reason behind satisfying the above mechanical properties (galilean in- 
variance, integrability, angular momentum, center-of-mass position) are understood to be 
the fact that the mass flux velocity Um is the same as the momentum density velocity Um- 
So we are verifying these properties with respect to the mass velocity rather than a volume 
velocity. All these demonstrations can be obviously reproduced taking the original set of 
volumed-based equations fll2ap to (]12dl) . 

E. Second law of thermodynamics and the apparent continuum limit volume/mass 
diffusion model 

For the second law of thermodynamics, momentum equation fl38bl) is rewritten: 

D 



2 



-f/^ ■ V ■ [pi + Hv] , (63) 



and then used to eliminate the macroscopic kinetic energy in the energy equation I38cl which 
becomes 

P-^ [e.J =t/,„ ■ V ■ bl + Hv] - V ■ [[pi + Hv] ■ U,] - V ■ {-KhVT) . (64) 
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Splitting the shear stress into its components due to Um and that due to J„, respectively, 
Hum and IIjv, equation flM|) can be rewritten: 



P-^ [Cin] =-pI: VUrn " ^Um ■ ^U^m " Hj, : Vf/^ (65) 

- V ■ [pi ■ J,] - V ■ [Hf/^ ■ JJ - V ■ [Hjv ■ J,] 

- V ■ (-K„VT) , 



where continuity equation fl38ap can be used to eliminate the first term on the right-hand- 
side. Then taking the definition of specific entropy s by the Gibbs local equilibrium equation 
equation (1^ becomes 

pT^ = - Uum : VUm - Hjv : VU^ (66) 

- V ■ [pi ■ J^] - V ■ [Uum ■ J^,] - V ■ [Hjv ■ J„] 

- V ■ (-K/.VT) , 



which may be rewritten as, 



or 



pT^ = - n„ : VUm (67) 

- V ■ [pi ■ J,] - V ■ [n, ■ J,] 



P^ + ^ V ■ [-K,,VT + i?T/€„Vp] = - : Vt/„ - ^ J, ■ V ■ (68) 



The last term on the right-hand-side of (1681) . as opposed to the first term, cannot be made 
to satisfy the positivity structure of the entropy production. Consequently, it appears that 
to fully satisfy positivity of the entropy production rate we must return to the full non- 
local-equilibrium structure of the volume-based hydrodynamic equations that involves the 
evolution of the additional microscopic structure evolution equation; this is the description 
presented in section [IIIBI We conclude that, while it may be difficult to justify a full consis- 
tency of a diffusive volume/mass model with the classical local equilibrium Gibbs equation, 
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beyond the local-equilibrium thermodynamics, there exists a volume/mass diffusion hydro- 
dynamic model satisfying all mechanical and thermodynamic properties. 
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